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We show that if the Rabi frequency is comparable to the Bohr frequency so that the 
rotating wave approximation is inappropriate, an extra oscillation is present with 
the Rabi oscillation. We discuss how the sensitivity of the degree of excitation to 
the phase of the field may pose severe constraints on precise rotations of quantum 
bits involving low-frequency transitions. We present a scheme for observing this 
effect in an atomic beam. 



It is well known that the amplitude of an atomic state is necessarily 
complex. The electric or magnetic field generated by an oscillator, on the other 
hand, is real, composed of the sum of two complex components. In describing 
semiclassically the atom- field interaction involving such a field, one often side- 
steps this difference by making the so-called rotating wave approximation 
(RWA), under which only one of the two complex components is kept, and 
the counter-rotating part is ignored. We discuss how the sensitivity of the 
degree of excitation to the phase of the field poses severe constraints on precise 
rotations of quantum bits involving low-frequency transitions. We also present 
a scheme for observing this effect in an atomic beam. 

We consider an ideal two- level system where a ground state |0) is coupled 
to an excited state and the <-> 1 transitions arc magnetic dipolar, with a 
transition frequency lo, and the magnetic field is of the form B — Bq cos{Lut + 
4>). We now summarize briefly the two- level dynamics without the RWA. In 
the dipole approximation, the Hamiltonian can be written as: 

H = eiao-ffz)/2 + git)a, (1) 

where g(t) = —go [exp(ja-'t + i(f)) + c.c] /2, ai are Pauli matrices, and e = u) 
corresponds to resonant excitation. The state vector is written as: 

\m) = (2) 



We perform a rotating wave transformation by operating on with the 

unitary operator Q, where: Q — (cto + crz)/2 + exp{iLut + i(f))(aQ — a^jl. The 

Schrodinger equation then takes the form (setting h. ~ 1): |^) = —iH{t)\£_{t)) 
where the effective Hamiltonian is given by: 

H = a{t)a++a*{t)a^, (3) 
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with a{t) = —{go/2) [cxp(—i2ujt — i2(p) + 1], and in the rotating frame the 
state veetor is: 

ic»)-Qiew) = (gg). (4) 

Now, by making the rotating wave approximation (RWA), corresponding to 
dropping the fast osciUating term in a{t) one ignores effects (such as the 
Bloch-Siegert shift ) of the. order of (go/uj), which can easily be observed in 
experiment if go is large EJu. Otherwise, by choosing go to be small enough, 
one can make the RWA for any value of uj. We explore here both regimes, 
and we find the general results without the RWA. 

From Eqs.^ and |[ we get two coupled differential equations: 

do{t) = iigo/2) [1 + exp(-i2wt - i2^)] Ci{t) (5) 

ti{t) = i{go/2) [I + exp(+i2tjt + i2(j))] Co{t). (6) 

Let |Co(0)p = 1 be the initial condition, and we proceed to find approximate 
analytical solutions of Eqs.5 and 6. Due to the periodic nature of the effective 
Hamiltonian, the general solution to Eqs.5 and 6 can be written in the form: 

E (J;)exp(n(-^2a;^-^20)). (7) 

Inserting Eq.^ in Eqs.5 and 6 we get for all n : 

a„ = i2nLoan + igo{K + ^^n-i)/2, (8) 
hn = i2nujbn + igo{an + a„+i)/2. (9) 

Here, the coupling between ao and bo is the conventional one when the 
RWA is made. The couplings to the nearest neighbors, a±i and b±i, are 
detuned by an amount 2uj, and so on. To the lowest order in (go/co), we 
can ignore terms with \n\ > 1, thus yielding a truncated set of six equations 
for do, &0j a±ij &±i- Now consider go to have a time-dependence of the form 
9o{t) = goM [1 — exp(— t/rsu,)], where the switching time constant Tsw is large 
compared to other characteristic time scales such as l/w and l/goM- Then, 
these equations can be solved by the method of adiabatic elimination, which 
is valid to first order in = (,go/4a;). (Note 77(i) — 770 [I — exp(— i/Tgu,)], 
where 770 = (goM/Aoj)). To solve the set of equations above, we con- 
sider first a_i and 6_i. Define /i± = (a„i ± 6_i), then one can write 
/i± = — i(2w + go/2)^± ± igoao/2. Adiabatic following then yields (again, 
to lowest order in i] ): a_i ~ and 6_i « 7700- Likewise, we can show that 
ai « —rjbo and 61 « 0. The only non-vanishing (to lowest order in rj with 
\Co{t = 0)p = 1) terms in the solution of Eqs.5 and 6 are: 

Coit) = co8ig'o{t)t/2) ~ 777e-*(2.*+2« sin(g^(i)i/2), (10) 
C,{t) = 7e-*('^*+«[sin(.9^,(t)V2) - 777e+'(2^*+2^) cos{g'o{t)t/2)], (11) 
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where g'Q{t) = 1/t go{t)dt = go [l - {t/Tg,,)-^ exp(-t/r„^,)] . 

To lowest order in ?y, this solution is normalized at all times. Note that if 
one wants to carry this excitation on an ensemble of atoms using a 7r/2 pulse 
and measure the population of the state |1} after the excitation terminates 
(at t = r when g' {t)t /2 = it /2 ), the output signal will be: 

|Ci(ff[,(T),0)|2 = ^ [1 + 2r?sin(2a;T + 2.^)] , (12) 

which contains information of both the amplitude and the phase of the driving 
field B. This clearly indicates that the Rabi transition probability depends 
on the total phase (fir = u)t + (j) oi the driving field. 
A physical realization of this result can be best appreciated by consider- 
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Figure 1. Left: Schematic illustration of an experimental arrangement for measuring the 
phase dependence of the population of the excited state 1 1): (a) The microwave field couples 

the ground state (10)) to the excited state A third level, |2), which can be coupled 

to |1) optically, is used to measure the population of |1) via fluorescence detection, (b) 
The microwave field is turned on adiabatically with a switching time-constant Tsw , and the 
fluorescence is monitored after a total interaction time of r. Right: Illustration of the Bloch- 
Siegcrt Oscillation (BSO): (a) The population of state as a function of the interaction 
time r , showing the BSO superimposed on the conventional Rabi oscillation, (b) The BSO 
oscillation (amplified scale) by itself, produced by subtracting the Rabi oscillation from the 
plot in (a), (c) The time-dependence of the Rabi frequency. Inset: BSO as a function of 
the absolute phase of the field with fixed r. 

ing an experimental arrangement of the type illustrated in Fig. 1 (left). And 
plot of the Rabi oscillations is shown in Fig. 1 (right). Under the RWA, the 
curve of Fig. 1(a) (right) represents the conventional Rabi oscillation. How- 
ever, we notice here some additional oscillation, which is magnified and 
shown separately in Fig. 1(b) (right), produced by subtracting the conventional 
Rabi oscillation {s\v?{g{t)/2)) from Fig. 1(a) (right). That is, Fig. 1(b) (right) 
corresponds to what we call the Bloch-Siegert Oscillation (BSO), given by 
r7sin((jrQ(r)T) sin(2(/)r). These analytical results agree very closely to the re- 
sults which are obtained via direct numerical integration of Eqs.5 and 6. Note 
that the BSO is at twice the frequency of the driving field, and its amplitude 
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is enveloped by a function that vanishes when all the atoms are in a single 
state. This osillation will be statLffer when the ratio of the resonance frequency 
to the Rabi frequency is large Wa u. One should keep track of the phase of the 
pyr'tptinn field at the location of the qubit Ofl for low energy qubit systems 
with fast driving. 

In conclusion, we have shown that when a two-level atomic system is 
driven by a strong periodic field, the Rabi oscillation is accompanied by an- 
other oscillation at twice the transition frequency. This extra oscillation can 
limit the qubit operations in Rabi flopping without a proper matching of the 
parameters. However, it has been shown that, this phase has potential applica- 
tion in teleportation. By making use of distant entanglement, this mechanism 
may enable teleportation of the phase of a field that is encoded in the atomic 
g^:^g)litude, and has potential applications to remote frequency locking 
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